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I. INTRODUCTION 



Integrable quantum field theories in d = 1 + 1 dimensions are more than just "theoretical 
laboratories". Apart from being very rich mathematical structures, they have proven to 
be of interest for a variety of problems in theoretical physics, ranging from string theory 
to statistical mechanics lattice models and problems in condensed matter theory At the 
technical level, the magic of integrable field theories can be traced back to the existence 
of an infinite number of non-trivial conserved charges M . These imply the conservation of 
individual momenta in many-particle collisions and the factorizability of the S-matrix. The 
two-body S-matrix satisfies the Yang-Baxter equation, in addition to the usual properties 
of unitarity, crossing symmetry and analyticity. In many cases, these properties of the S'- 
matrix are so restrictive that one may set up a bootstrap program and determine the exact 
S- matrix, up to the famous "CDD ambiguity" ||. 

One special class of integrable field theories, that we will be studying in this paper, are 
the integrable massive perturbations of conformal field theories. Such theories are obtained 
by starting with a conformal field theory (CFT) and perturbing it by a relevant operator that 
preserves integrability. We will see the technicalities later. The ultraviolet (UV) properties 
of the CFT will not be affected by the perturbation, but the infrared (IR) behavior will 
change. The modified IR behavior can be conformal or massive (finite correlation length). 
In the latter case, the IR behavior is given by a massive particle theory with f actor izable 
scattering and one may employ the above-mentioned bootstrap program to identify the exact 
S'-matrix. 

Once the S-matrix of an integrable field theory is known, one may study its thermody- 
namics by following a procedure called Thermodynamic Bethe Ansatz (TBA) [|,[|. In the 
context of applications in condensed matter problems (such as the scattering of edge cur- 
rents in the fractional quantum Hall effect, ||), the TBA results offer concrete predictions 
that can be tested in experiments or simulations. Another nice feature is that, among the 
quantities that can be computed by the TBA are the central charge and the scaling dimen- 
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sions of the UV limit of the theory. In this way, the TBA may be used as a check on the 
validity of a conjectured S'-matrix for a perturbed CFT. The name of the TBA technique, 
which will be reviewed briefly in section III, is a little misleading, since the TBA is an exact 
non-perturbative procedure and not a trial-and-error method. 

The theories that we study in this paper are examples of integrable field theories with 
super symmetry. Adding supersymmetry is of course natural in the context of string theory 
applications, but even in the context of statistical mechanics lattice models supersymmetry 
is meaningful. For example, the tricritical Ising model in two dimensions 0] is an example of 
a physical system that realizes superconformal symmetry. If we start from a superconformal 
field theory and choose a perturbation that preserves supersymmetry we end up with a 
massive theory with supersymmetry. In such a case supersymmetry may be added to the 
bootstrap ingredients, and one may try to identify exact S-matrices by exploiting both 
supersymmetry and factorizability. 

In this paper we will study the TBA for a series |3J] of N = 1 supersymmetric scattering 
theories with 2n particles, n some positive integer, arranged in n supermultiplets of one 
bosonic and one fermionic particle each, with masses labeled by an integer k and given by 



where j3 = l/(2n + 1) and k — 1, 2, • • • , n. It has been conjectured |4j] that these scattering 
theories correspond to an integrable perturbation of specific superconformal field theories of 
central charge c n = — 3n(4n + 3)/ (2n + 2). Note that this is a series of non- unitary theories. 

One important observation is that the S'-matrix for these particular theories is non- 
diagonal, since we can have, for example, a scattering process of the form bibj — > fifj, where 
bi{fi) is a boson (fermion) in the z-th supermultiplet. This makes the TBA analysis more 
difficult than in the case of diagonal scattering, where the thermodynamic limit can be 
studied via a set of coupled integro-differential equations. In the non-diagonal case we will 
have to diagonalize a tranfer matrix in order to obtain a tractable system of equations. 

In general, the transfer matrix that features in the TBA for non-diagonal scattering 



sin(fc/37r) 
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theories is identical to the transfer matrix of a corresponding statistical mechanics lattice 
model, which is obtained by interpreting the two-body S-matrix elements as Boltzmann 
weights. In the case of the above N — l supersymmetric theories, this lattice model is of eight- 
vertex type, and one might expect that the analysis would be quite involved. However, it has 
been observed that in a rather general setting || the eight-vertex models that correspond 
to N = 1 scattering matrices satistfy a technical condition callled the free fermion condition. 
This case is no exception, and we shall use this special property to determine the eigenvalues 
of the transfer matrix and to complete the TBA analysis. 

One of the motivations for us for deriving our N = 1 supersymmetric TBA systems has 
been a proposal made by E. Melzer in ||. He conjectured that the N = 1 TBA systems 
can be derived by starting from scattering theories with N = 2 supersymmetry, and twice 
as many particles, and applying a procedure called "folding" on these N = 2 TBA system. 
These N = 2 TBA systems were first studied by P. Fendley and K. Intriligator in JFj. Melzer 
checked that the "folded" TBA systems had all the right properties to be related to the 
N=l superconformal theories that underly the N=l S-matrices, but was unable to give a 
direct derivation starting from the N = l S-matrices. In this paper we fill in this missing link 
and confirm Melzer's conjecture by showing that the true N = 1 TBA systems are indeed 
identical to the folded N = 2 systems of ||. 

In earlier work, the bosonic scattering theories that underly the N—l and N = 2 
supersymmetric scattering theories that we just discussed have been related by a similar 
folding, both at the level of the TBA and at the level of the S'-matrices |8]|| . Our results in 
this paper clearly suggest that it should be possible to find a "folding relation" between the 
N = 1 and N = 2 supersymmetric theories, directly at the level of the S'-matrices. We will 
address this issue in a forthcoming publication [ |TTf . 

This paper is organized as follows. In section II we review the formalism for supersym- 
metric particle theories with factorizable scattering and we introduce a specific series of 
N = 1 supersymmetric theories. In section III we discuss the TBA technique and highlight 
some aspects that will be relevant later. In sections IV and V we obtain the TBA system 
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for the N =1 supersymmetric theories and verify the UV and IR limits. In section VI we 
present Melzer's folded N = 2 systems and demonstrate that they agree with our N — 1 
systems. Our conclusions are formulated in section VII. 



II. SUPERSYMMETRIC PARTICLE THEORIES IN 1+1 DIMENSIONS WITH 

FACTORIZABLE SCATTERING 

In this section we briefly introduce our n-component supersymmetric scattering theories 
and establish some notation. For a more detailed discussion, see j4|. 

We are interested in perturbed superconformal field theories decribed by an action of the 
form 

S X = S + \J (GUGLi<M d 2 x , (2.1) 

where S is the action of the unperturbed theory, which we take to be the minimal supercon- 
formal model -M(2, An + 4) of central charge c n = — 3n(4n + 3)/(2n + 2). The perturbation 
contains 4>h,hi which is a primary field in the Neveu-Schwarz sector of the left and right-chiral 
superconformal algebras. The perturbations in (|2.1|) are manifestly supersymmetric. If we 
consider the case with h = /i(i,3) we obtain an integrable theory. This is what we meant by 
massive deformations of CFT's in the introduction. It can be shown by means of a slight 



variation of Zamolodchikov's counting argument |12[ that for this perturbation there is at 
least one non-trivial integral of motion, which ensures the factorizability of the S'-matrix. 
Since the perturbation has the dimension of action, we are automatically introducing a scale 
in our model and the resulting theory is not conformal invariant. On the other hand, in the 
ultra-violet (UV) limit we obtain a (super) conformal theory. This UV limit is simply the 
theory described by S (the unperturbed theory) since in the limit of extremely high energy 
scattering, the particles do not "see" any finite energy scale. The infra-red (IR) regime of 
the theory is completely described by a factorizable scattering theory, which, according to 
the conjecture in H can be described as follows. 
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The IR scattering theory contains 2n massive particles, arranged in n supermultiplets 



(bi,fi), with masses m ; as given in ( |1.1| ). We denote by Aj(^) be any particle, boson 
or fermion, from any multiplet, with asymptotic momentum p® = mjcosh(0j) and p\ = 
m;sinh(#j). Multi-particle states are then written as 



I -^-h ) ^-il (^12 ) • • -^ijv (^ijv ) ) in(out) 



(2.2) 



where 9^ > 9i 2 > . . . > 9 iN for in-states and the other way around for out-states. 

The two-body S-matrix will carry two labels, % and j, that tell us which supermultiplets 
the particles in the scattering belong to. From Lorentz invariance it is easy to see that the 
two-body S-matrix can depend only on the rapidity difference 9 = 9^ — 9j. The complete 
two-body S-matrix is written as a product of two pieces 



S^{9) = SH{9) S™(0) 



Mi 



(2.3) 



where S^p is a piece of the S'-matrix that mixes bosons and fermions and Sffl is a bosonic 
S'-matrix that acts only in the bosonic sector. 

For the model at hand, the bosonic piece Sffl is the S-matrix found by P. Freund, 
T. Klassen and E. Melzer in [T3j and given by 



Sf{9) = -F^ p {9) F { ^\ +m {9)...F { 



1 2 



(9) F {l+J)p (9) 



(2.4) 



with F a (9) = ^Si-isin& In a basis § iven b y IW' l^/i)' l/A">> \f*fj) the S bf P iece is 



of the form 



sinh(S)— i sin(a7r) 



S f F (9) = fW{9) 



\ 



^ 1 ^ 



+ g m (0) 



, (2.5) 



1-tt -i(t + t) 

-t + i 1 + it 

l + it t-i 

, -i(t + t) 1-tt j yuuu-iy 

where t = tanh((0 + log(mj/mj))/4) and t = tanh((0 — log (mi/my ))/4). The functions 
/M(0) and c^'](0) are related by 



10 
10 
-1 
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/M(0) = ^^m. 



a 



Ai 



2 cosh(fl/2) + (p 2 + p- 



(2.6) 



cosh(#/2) sinh(#/2) 

where p = (rrii/mj) 1 ^' 1 and a = — sin(/37r). This specific form of the Bose-Fermi S'-matrix 
Sgp is almost completely fixed by N = 1 supersymmetry (which allows only two free functions 
/™(0) and g™{0)) and by the Yang-Baxter equation, which fixes the ratio of /™(0) and 
<7^(0) up to one free constant a. The value of a is then dictated by consistency with the 
bound-state structure induced by the bosonic factor Sffl in (|2.4j ). 

The function g^(9) is constrained by the conditions of analyticity, crossing symmetry 
and unitarity, 



analyticity : 



crossing symmetry : 



unitarity : 
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9 [ij] *(0), 



g&\m-6) = gM(6) 



(2.7) 



gM(9)g 



1 , 2 /sinh 2 (e/2)+i(p 2 +p- 2 ) 2 

1 + a 2 m imj f cosh . (e/2) 4 sinh . (e/2) 



-i 



An explicit integral expression for g™(0), which will be important for the TBA analysis, 
can be found as follows. If we define the angles Si = \{i + j)/37r and 62 = — (i — j)/?7r) 
we can write the unitarity equation as 
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M(a\ n M 



sin(f -f)sin(f + f)sin 2 (f; 



sin(f - Si) sin(f + Si) sin(f - S 2 ) sin(f + <f 2 ) 



(2.8) 



The solution for the one component case was found in [14 where the equation for g{6) is: 

sinh 2 ^) 



9A(e)g A (-e) 



sinh 2 (|) + sin 2 (A7r) 



(2.9) 



where A is some parameter. The solution is then 

sinh(f) 
smh(|J + 1 sm(A7r) 

We can now find a solution for 



00 Gftsinh(At)sinh((l - A)t) ,tB 
exp I % 1 - — Arm ... sm( — 
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(2.10) 



t cosh 2 (|) cosh(t) 
by writing a product of three functions of the form 



gM{9) = 9aM9aM 

9A-A 6 ) 

where A 1 = \(i + j)f3, A 2 = |(1 - (i - j)0) and A 3 = \. This is the form of g [ij] (0) that we 
shall use in later sections. 

Before closing this section, we remark that the S"-matrix S^p in (|2.5| ) is related to the 
hatted matrix S^ F of reference by Sbf = nS"^, where II is an ordinary (not graded) 
permutation. Making this choice (rather than working with IT gra ded) effectively compensates 
for minus signs that would be induced by the fermionic nature of the particles fc, and it will 
allow us to treat all particles as bosons when performing the TBA. Note however that this 
choice breaks the manifest supersymmetry of our formalism. 



III. SOME ASPECTS OF THE TBA 

In the previous section we have recalled the conjecture, made by one of us some five 
years ago, that the supersymmetric and factorizable S-matrices (|2.3| )- ([2.5| ) are actually the 
true S-matrices for the perturbed conformal field theories ( |2.1| ). This claim was made on 
the basis of the compatibility of the set of conserved quantities (which can be derived in the 
perturbed CFT framework) and the bound-state structure that is implied by the proposed 
^-matrix. One way to obtain further support for the conjectured identification is by using 
the TBA. In brief, the TBA makes it possible to express some of the UV CFT data (such 
as the central charge) in terms of the 5*-matrix data, providing a non-trivial check on the 
correctness of the proposed S'-matrix. Before we come to this analysis, we shall in this 
section briefly review the TBA for diagonal and non-diagonal S'-matrices. 

If we want to know the energy levels of a free particle on a circle of radius R, we just have 
to send the particle on a round trip around the circle and impose a matching condition for 
the wave function. This will constrain the possible momenta and from that we can compute 
the energy levels. In an equation this reads simply 

e ip2wR = 1 . (3.1) 
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From here we obtain p n = n/R for some integer n and energy levels given by E n = p 2 l /2m = 
n 2 /2mR 2 . 

Integrable models are not free of course, but due to the fact that multi- particle scattering 
can be studied two particles at a time, and that this scattering will be purely elastic, we 
obtain enormous simplifications and the parallel with a free particle on a circle is very 
appropriate. The only input we need is the exact ^-matrix and the mass spectrum of our 
theory. Let us explain the case where the S'-matrix is diagonal first, just to make the concepts 
more clear. A more complete discussion of these matters can be found in ||U . 

3.1 The diagonal case 

The idea of the TBA is very simple physically. We consider a situation where M particles 
are put on a circle of lenght L = 2ttR. We then take one particle, of mass m,i and rapidity 
6i, and send it on a round trip around the circle. Since we the S'-matrix is diagonal we will 
simply pick up a phase every time our particle meets one of the remaining M — 1 particles. 
The wave function should come back to its original value after the particle has completed 
its trip and this implies 

e im ' sinh(e ' )L J] S [lj \e tJ ) = 1 . (3.2) 

Taking the logarithm of this equation and going to the thermodynamic limit (M — > oo, 
L — > oo, with M/L constant) we obtain the following equation 

27rP;(0) = m;Lcosh(#) + J2 J d6 ' Pj( e ')&j( e ~ e j) > ( 3 - 3 ) 

j 

where P%{9) is the density of available rapidities for particles of type i, and Pi(6) is number 
of such rapidities that are actually occupied. The phase shift 4>ij(0) is given by 

<t>M = -i qq^ 2 - ■ (3-4) 

Let us define the pseudo-energy levels €i(9) by 



The Pi(9) are such as to minimize the free energy. From this, we get an equation for the 

/rift' 
—^(e-e^Hi + e-^) , (3.6) 

The total energy of our system can be written in terms of e,(#) according to 

E(R) = - V ^ [dft cosh(^) ln(l + e-^ {e) ) , (3.7) 

i 27T J 



We shall call the equations (|3.6|) , (|3.7| ) the TBA equations of the scattering theory. 

From the TBA equations we can find numerically the Casimir energy for any R. In the 
UV limit we can do even better and use the amazing Roger's dilogarithm identities to derive 
closed- form expressions (an explicit formula will be given in the next subsection). 

3.2 The non-diagonal case 

In this case we have conceptually the same situation. The difference is that now a 
particle, when taking a round trip, scatters non-diagonally with all the others so that the 
periodic boundary conditions lead to an equation involving a transfer matrix of size 2 M x2 M . 
This transfer matrix can be written as 

( 7 jiu,.,-.u)y^-^ m ft 2 ...ft M) = 
£(,s^)#(0 _ e 1 )(s^y r f s i(e -e 2 )... (s^y^_ lSM (e - e M ) , (3.8) 

In this notation, the indices j and ji denote the different multiplets = 1,2, ... ,n) and 
the indices ri, si and ti denote the two states b and / within each multiplet. The periodic 
boundary condition is now expressed as (see, for example, 0) 



3 imj ( sinh(#;)L 



£ (Tin\nn-m)^ ^ = _^ ? / = 1, 2, . . . , M , (3.9) 



where ip is a 2 M -component wave-function. 

To make further use of this equation, we have to diagonalize the transfer matrix, or 
at least find some relationship among the eigenvalues that will allow us to write the TBA 
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equations. In general, the problem of diagonalizing a transfer matrix is extremely difficult 
and many techniques have been developed to handle such problems O]. We will, in the 



next section, use a method developed by B. Felderhof [T6|] to diagonalize the transfer matrix 
for the eight-vertex model in the case where the Boltzmann weights satisfy the free fermion 
condition. Remarkably, the resulting TBA equations can be cast in a form identical to ( ^.6| ), 
Q3.7|) , this time for a set of particles that the contains the original set plus an additional 
particle of mass zero (see section V). This shows that the form ( |3.6|) , ( |3.7|) of the TBA 
equations is universal, applying both to diagonal and non-diagonal cases. 

Returning to the Casimir energy, it can be shown that the UV limit is given by 



E = E{mR -> 0) 



— v 



X ; 



a 



Vi 



l + Xi 1 + yi _ 

where C(x) is the Roger's dilogarithm function, which is defined by 

1 r x 



(3.10) 



C(x) 



2 Jo 



dt 



\n(t) ln(l - 1) 



t 



t 



(3.11) 



and Xi = e and y% — e The yi are the solutions of 



X i 



0i'= IK 1 + 



(3.12) 



where the Nij are given by 



i -J 



2tt 



(3.13) 



The primed indices %' refer to the massless particles only. 



IV. DIAGONALIZING THE TRANSFER MATRIX 

As we have seen in the previous section, we need to gain control over the eigenvalues of 
the trace of the transfer matrix T(9\9x, 9n)- Notice that this is a inhomogeneous transfer 
matrix, depending on all the differences (8 — 9{). Due to the Yang-Baxter equation, tr[T(#)] 
and tr[T(0')] commute for any 9 and 9' . This means that the eigenvectors of tr[T(0)] are 
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^-independent. The eigenvalues will of course depend on 6 and this dependence is precisely 
what we need to know to do the TBA. 

In section II, we mentioned that the non-diagonal part of the S'-matrix, called S^p, 
can be derived by imposing supersymmetry and the Yang-Baxter equation. A remarkable 
property is that the elements of Spj* F , when interpreted as the Boltzmann weights of a lattice 
statistical mechanics model, satisfy the free fermion condition. Explicitly this means that, 
for sf] given by 



S 



m 

BF 



a + 








d 





b+ 


c 








c 


6_ 





d 








a_ 



(4.1) 



the a + , a_, b + , 6_, c and d satisfy 

a + a_ + = c 2 + d 2 



(4.2) 



It can be proved that an eight-vertex model that satisfies this condition can be written as 
an XY model with a magnetic field. The magnetic field is given in terns of the Boltzmann 
weights 



H 



a 2 _ + bl-a 2 ,- b 2 _ 



(4.3) 



2(a + b- + a_6+) 

Using the results from section II we find that H = — 1 for all i, j. This is the critical point 
of the XY model. 

For the determination of the eigenvalues of tr[T(0)] we will follow a method given by 



Felderhof pq] . We shall explain the idea of this method and apply it to our case. The 
strategy will be to find an inversion relation for tr[T(6 1 )] which can then be used to write 
a functional equation for the eigenvalues. It will turn out that, rather miraculously, the 
analysis for general n can be done by extending the analysis for n — 1, which was presented 
in |T7||. Our presentation in this section therefore follows [T7| rather closely. 

It is not difficult to see that the S'-matrix Q2.5|) satisfies the following relation 
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s 



M t 

BF \ 



+ in) 



VK 



gM(6) 



a_ c 

-6_ -d 

-d -b + 

c a. 



(4.4) 



where we are using a obvious notation. If we make a similarity transformation by means of 
a matrix M (acting on the indices r = b, /)[] given by 

1 



M 



(4.5) 



we obtain the following 



gM(6 + 



in) 



gM(B) 











-d 





a + 


c 








c 


a_ 





-d 








-6_ 



(4.6) 



/ 



It is easy to S66 t licit S again satisfies the free fermion condition. 

>BF 



Returning to the transfer matrix (|3.8p , let us write Si = S 



9i), etc., so that 



tr[T(0)] = tT[SxS 2 ...SM} and similarly tr[T(0 + nr)] = tr[T(0)] = tT[3 1 3 3 ...3 M ]- The product 
tr[T(0)] tr[T(0 + in)] can then be written as 



tr[T(0)] tr[T(0 + in)] = tr (^ ® <g> 5 a ) . . . (SW <g> Sat) 



(4.7) 



where the notation ® denotes a matrix product of the s-indices and a tensor product on 
the r-indices. It can be proved that there is a similarity transformation with a constant 
(^-independent) 4x4 matrix X (acting on the indices r h fj) that puts Si ® Si in a triangular 



1 Notice that the indices r and s on the two-body 5-matrix (5*^)^/ are not on equal footing in 
the transfer matrix fl3,8|) . In the usual terminology, the r are called matrix indices and the s are 
quantum indices. 
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form. We do not give the explicit form of X here, but mention the important fact that it 
depends on j, ji through a parameter 0, given by 

2cd 



tanh(0) 



(4.8) 



a + b + + a_5_ 

In terms of the elements of S^p we get tanh(0) = —arnj, i.e, it only depends on the index 
j, which is common to all factors in the transfer matrices. We can thus use one and the 
same X to put all the Si <S> Si in triangular form, and the total trace tr[T(0)] tr[T(8 + is 
invariant under this operation. We find 



-sf$(p-e l )®sif(e-o l 



Mix. 



M + 


* 


* 


* 







* 


* 










* 











M_ 



/ 



where M + , M_, F + and F_ are given by 

M + = a + a_ — d 2 , M_ = a + a_ — c 2 
F± — ± ^sinh 2 (0)a_|-6-|- + cosh 2 (0)a T 6 T =)= 2sinh(</>) cosh(0)c<i 

which can be rewritten as 
-4# 2 



M+(0) 
M_(0) 
F + (6) 
F-(0) 



sinh ( 
-4^2 ( 



sinh 2 ( 
-4^2 ( 



sinh 2 ( 
-4^2 ( 



sinh 



cosh( ) cosh( ) smh( ) smh( , 

. 1 . 6 1 + 7r . 8-i/3-ir. .8 + i/3 + TT. .6 - if3 + n , 
sinn( ) smh( ) cosh( ) cosn( 



. 6-if3.Tr 
sinn( J cosn( 



2 ' v 2 

8 + z/3_7r N , ,0 — if3 + n , 
) smh( 



cosh( 



+ i/3+7i\ 



cosn( ) smh( ) cosn( ) sinn(- 



(4.9) 



(4.10) 



(4.11) 



act 



2 , v 2 v 2 y v 2 ' 

where /3± = (jz ± j)/3 and #2(0) = (0) (0 + in). The entries of the matrix 
diagonally on the indices si, with the ± sign equal to +1 (—1) if total in-state is bosonic 
(fermionic). 
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Given this triangular form for Si ® Si, it is very easy to compute the eigenvalues of the 
product tr[T(#)] tr[T(# + iir)]. Denoting the eigenvalues of tr [T(9)] by A(9), we find 



K{9)A{9 + m) 



-1 



M 



M 



11 M, 



\i=i 



9 l ) + l[M49-9 l ) + F 

1=1 



■ M M 

]jF + (9-9i) + l[F-(9-9 l ) 

1=1 i=i 



(4.12) 



where F is +1 for a bosonic state and —1 for a fermionic eigenstate. Notice that the 
dependence of M± and F± on j, ji has been suppressed in the notation. In terms of reduced 
eigenvalues X(9), defined by 

Oi) 



M „[jh)( 

the equation (|4.12j) can be rewritten in the following factorized form 



(4.13) 



X(9)X(9 + in) 



F 11 11 cosh( )smh( 1 

k=n h =i A z 

. f[ sinh(^^^)cosh(^±^±l^)' 
k=llk=l 2 2 

pfrir ^ e-e lk +i{k-j)fr ^ . e-e lk -i{k + j)pn. 

b || || cosh( - ) smh(- 

k=i i k =i 



) 



(4.14) 



2 ' 2 

■ Air • <h.~M-j ) 
+ || || smh( ) cosh( 

k=l l k — ~" 

where we assumed that out of the M particles on our circle, Mj, are of type k, J2k ^k — M. 
Recall that the auxiliary particle that makes the roundtrip is of type j. 

The X(9) are 2iri periodic meromorphic functions, and as such they are completely de- 
termined by their zeros and poles. Let us denote by z + the zeros of the first factor in the 

rhs of (|4.14j ) and by z~ those of the second factor. They satisfy the equations 

>z--e lh +i(k~j)i3TT, 



tanh( : 



nn , 

k=n k =i tanh( 



z+-e t +i(k+j)fS-n 



-F 



ft ff taUl ' ( 



k=i i k =i 



Denoting by x m the real solutions of the equation 

11 11 . ,x—B u +ikf3-K, ~ 

k=n k =i tanh( \ ) 



z --e lk -i(k+j)pn 



tanh( 



-F . 



(4.15) 



(4.16) 
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one may check that the following satisfy ( |4.15| ) 

1 — e + _ 1 — e~ 

Z t+,m = X ™ ~ ijpTT + Z 7T , Z £ _ )m = X m + ij ^7T + i 71, (4.17) 

Out of each pair each of zero's zf ± , e 1 * 1 = —1, 1, one will come from \{9) and the other 
from \(9 + 7rz), which gives us two choices. It turns out that the choices to be made for the 
two factors (i.e. the choices for e + and e~) are correlated (see, for example, WA), so that for 
each m there are two possible factors contributing to \{6) 

sinh(^|^) sinh(^|^) , e = -l,l. (4.18) 

We thus arrive at the following expression for the eigenvalues 

^{x m ,e m }(0) = Constant x JJsinh( - * m + iej^-) cosh( ^ * m - iej^-) , (4.19) 

The counting of the eigenvalues (there should be 2 M of them), works out as follows. 
There are M — 1 solutions x m in the sector F — — 1, which, with freely chosen e m , gives 
the correct number of 2 M_1 eigenvalues in that sector. In the sector F — 1 there are M 
solutions x m , and in order to obtain the correct number of 2 M_1 eigenvalues there, we need 
one constraint on the allowed values of e m . Explicit results for some small values of M 
suggest that this constraint takes the form T[ m e m = 1. 

The expressions (|4.13|) and ( 4.19|) for the eigenvalues A ((9) will be important ingredients 
in the TBA analysis, which we present in the next section. 



V. THE THERMODYNAMIC BETHE ANSATZ 

5.1 The TBA system 

In order to do the thermodynamics we have a large number M of particles on a circle 
of lenght L. We send one of the particles, of type j and rapidity 9, on a round trip and 
scatter it off as it moves, one particle at a time due to integrability. When it comes back to 
the original point, we impose a periodic boundary condition. From this condition we may 
derive the allowed rapidity configurations, and from there the thermodynamics. 
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As we have seen before, the full S- matrix can be written as a scalar factor S^(6), given 
by (Ell), times the 4x4 matrix S l ^(9). We will denote the product S% j \9)g^\9) by Z^\6). 



The periodic boundary condition ( |3.9|) then takes the following form 

We denote by pk{9) the density of occupied states and by Pk{9) the density of available 
states for particles of type k, k = 1,2, ...,n. Similarly, we introduce the symbol Pq{6) 
for the density distribution of solutions x m to the equation ( |4.16|) . This is further split as 



Pq{9) = Pq{6) + Pq{6) where the densities p , p refer to those x m for which the corresponding 
e m = +1,-1, respectively. Notice that po(9) can be viewed as a density distribution for a 
new type of "particle" with label 0. This explains that the final TBA system will be in 
terms of n + 1 rather than n types of particles. 

From (|4.16 ) we obtain the density distribution of the available rapidities for the "parti- 



cles" of type 

, , " r , , , ls d , ftaDh(U9-0 , -ik/3Tr))\ , . 

fcrftW = g/^W^ln {-^__AJ^ . (5.2) 

Notice that the "mass term", which in general occurs on the r.h.s. of TBA equations of 
this sort, is absent. This means that the auxiliary particles of type should be viewed as 
massless. 

Taking the derivative of the imaginary part of the logarithm of (|5.1|) leads to 



271-^(0) = TTij cosh(0)L + J dd' 



X-( l i^ k \9-9') 



k=\ 



89 lsinh(0-0') 



+p o (0')^Imln(A£(0 _ e>)) + p o (0')^Imln(P o (0 - 0')) 



(5.3) 



where Aq, Aq are defined by 



A J = sinh(i(0 + ijfiir)) cosh(i(0 - ij/3ir)) 

\i = sinh(i(0 - ij/3*)) cosh(^(0 + ij/3*)) , (5.4) 
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Let us define some quantities that will make the notation lighter. Define Q% k \9) and 
&(6) as 



' w m Uinh(9) I 



&(0) = 2— Imln(A^)) = -2— Imln(A^)) . (5.5) 

Notice that the same quantity occurs in the r.h.s. of the equation ( |5.2| ). 

Using Pq{6) = po(9) + Pq{6) we can now eliminate po(9) from fl5.2|) , ( |5.3|) to obtain 



2nPj(6) = m jC osh(9)L 

+ E /^ / (p fe (^($fk^-^-^(^^^)(^-^)+Po(^^^-^) , (5-6) 
where we introduced the convolution 

(^'*$ fc )(0) = — / d#'^(#-#')$ fc (#') • (5-7) 

27T J-oo 

We can now vary the Pi(6), Po{9) so as to minimize the free energy. In terms of the 
quasi-particle energies e a (9), a = 0, 1, 2, . . . , n, defined as 

Pa{9) e~ e ^ 

PM ~ l + e--W ' (5 - 8) 

this gives the equations 

6j(9) = mj Lcosh(9) - ]T l(®z k] ~ -($ j *$ fc ))*hi(l + e - £fc )J (9) 

k=i ^ ^ ' 

+ ($ J '*ln(l + e~ e °))(#) 

n 

e o (0) = -£(<& fe *hi(l + e - e *))(0) , (5-9) 

k=l 

which have the general form ( |3.6|) . The total energy is given by (|3.7|) with the sum running 
over the index a = 0, 1, 2, . . . , n. 

Now that we have the TBA equations we can explore the UV limit and obtain analytical 
solutions for the ground state energy. This will be done in the next subsection. 



5.2 The UV limit 

The UV limit provides an important test for the conjectured 5-matrix since we can 
compute the ground state energy analytically in this limit and compare it with the value 
predicted by CFT: E = — ^c e ff, with c e g = c — 12(/t min + h m m), where h min (h m i n ) is the 
smallest conformal dimension in the theory. In the case of unitary theories h m i n = h m i n = 0, 
corresponding to the identity operator, which is always present in the operator content of a 
unitary CFT. In the case of non-unitary theories we have operators with negative dimensions 
and the effective central charge picks up the term — 12(/i min + h m i n ). 

In the CFT's that we started from in section II we have c = c n = —3n(4n + 3)/(2n + 2) 
and h min = /2, min = — | so that c c s = 3n/ (2n + 2). 

To obtain the ground state energy in the UV limit from the TBA developed in the pre- 
vious section we shall work out the general result ( |3.10|) - (|3.13|) for these particular theories. 

We denote by Ajk and Bk the following basic integrals 

/■oo Afj poo Afj 

A»= ^ *?(*). *i = ^{0). (5.10) 



4 J* - / TT! *i \P) > D i 

J —oo Z7T j — oo 

Using the explicit form of & z k {6) it is easy to compute Ajk- We have 



sinh(#) J 



= 1*4 In ( ^ {e), ^ H ^{ 6) ' Sinh( ^ )] e * W> <ffl { Q)\ , (5.11) 
69 \ [gA 3 {0)/ smh{9)\ j 

where F(t) is a smooth function that falls exponentially to at t — > ±oo and #a is the 
prefactor of formula (|2.10|) . From equation (|2.11| ) we see that the u g&" piece gives ~ + | — ~. 
The exponential gives zero since after the derivative with respect to 9 we get an integral of 
a smooth function F(t) times sin(^) and by the Riemann-Lebesgue lemma this will be zero 
in the limit 9 — > ±oo. Each factor F a (with a^O) from the bosonic S"-matrix ( |2.4| ) can be 
easily seen to give —1, and therefore we obtain 

/°° d9 f) r-n 
--l m \n(S [ i k] (9)) = -2k + 8 j>k , j > k . (5.12) 
-oo Z7T Ou 
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Collecting these results we see that Ajk = —2j + 8j t k + \ for j > k. The computation of Bj 
is straightforward and gives B 3 ■ = 1. This is all we need to write the equations ( |3.12|) for the 
pseudo-energies in the UV limit 

n n 

x 3 = {l + x ) N °>* l[(l + x k ) N - k , x = l[{l + x h ) N °* , 2/0 = 1- (5.13) 

k=l k=l 

with Xj = e~ ej (°\ x = e~ to( -°\ y 3 = e~ tj ^ and y = e~ t0 ^°°\ Note that due to the presence 
of the mass term, in the limit 9 — >• oo the €j(0) diverge and this implies directly that y 3 = 0. 
The N a> b are given by Nj t k = Ajk — \B 3 B k and N j = Bj so that 

N j)k = -2k + 5 hk , N d = 1 , iVo,o = . (5.14) 

for j > k. The other N a ^ are obtained using the symmetry N a ^ = N b<a . 
The (far from trivial) solutions for equation (|5.13|) are 



sin 2 ! 1 ^) sinf 37r 



sm ( L fcr) sin ( i i+f ) sm («) 

We can now use the Roger's dilogarithm technology and compute, using equation fl3.10|) , 



the effective central charge c e fj 



6 



C e ff = 



4n+4 



, , :m 2/ (2j+l)7r x I I 4 cos 2f 7r 1 / I 9 

3=1 V Sin (_ 4n 4 J / \ ^ cub V4n+4^/ V " 



3n 

(5.16) 



2n + 2 ' 



in agreement with the value in the original conformal field theory. 

We found the solutions Q5.15p for Xj, x and the dilogarithm identity (|5.16p in the second 



paper of J7J, which deals with TBA systems for N = 2 supersymmetric theories. The relation 
between the N = l and N = 2 TBA systems is the topic of our next section. 



VI. MELZER'S FOLDING 

In a very beautiful paper , E. Melzer studied a series of new TBA systems with some 
mathematical applications in mind, namely, the study of number theoretical identities of 
the Gordon-Andrews type. He developed some TBA systems that provide such identities. 
Melzer's obtained his "first TBA system" by a folding in half a TBA sytem derived by 
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P. Fendley and K. Intriligator |7j for specific N = 2 supersymmetric theories. He then 
conjectured that this folded TBA system should correspond to the N — 1 supersymmetric 
scattering theories of which we discussed in section II. Now that we have derived explicitly 
the N — l TBA systems we are in a position to check, and confirm, Melzer's proposal. 
Let us briefly summarize Melzer's discussion. We consider a general TBA system 

m n r +oc d9 



E(R) = - V ^ / °° ^ cosh(tf) ln(l + e-^W) (6.1) 

" 2,71 J-oo 2lX 



with 



e a {6)=m a Rcosh{8)-'£(<S> a>b *ln{l + e- e >>))(6) , (6.2) 

b 

where the kernels Q a ,b are symmetric in a, h. Let us now assume that we have 2n types of 
particles, a = 1, 2, . . . , 2n and that the TBA data m a , $ a fe possess the following symmetry 

m a = m 2n +l-a , $a,6 = $2n+l-a,2n+l-b • (6.3) 

In such a situation we can define a folded TBA system for half the number of particles, 
a = 1,2, ... ,n, by defining a folded kernel 

$5 ded = * a,b + $a,2n+l^ , O, 6 = 1, 2, . . . , 71 . (6.4) 

The N = 2 supersymmetric TBA systems of contain 2k — 2 massive particles and 
auxiliary massless species labeled as 0, 0. They correspond to massive scattering theories 
that are obtained by starting from the N = 2 minimal superconformal field theories and 
perturbing by the most relevant supersymmetry preserving operator. This corresponds to 
the Landau- Ginzburg superpotential 4^ — ^ = 2, 3, .... The masses are given by 



m *= Z7,*-\ » a = 1,2,..., 2* -2 , (6.5) 



sin 



2fc-l/ 

.N=2 fl,„ faN=2 



so that m a = rri2k-i- a , and the kernel has the symmetry property 

2A-1-6 ( we identify the labels = 2A; — 1). We can thus use ( |6.4j ) to define a 



folded version of the N = 2 TBA system. According to Melzer's conjecture, this folded N = 2 
system will be the same as the N—l system that we derived in this paper. 
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Identifying n = k — 1, we see that the mass spectra ( |1.1|) and ( |6.5|) indeed agree. We 
also checked that 

< b =1 (d) = K b =2 (0) + <2n\i- b (d) , a, b = 0, 1, 2, . . . , n , (6.6) 

where $^ 2 (^) are the kernels obtained in J7| and are the kernels in our TBA 

system ( |5.9|) . For the kernels $fc,o(#) this result is elementary while for the kernels $^^(9) 
some tedious cosine Fourier transforms are needed. 

In the UV limit the folding relation can be expressed as 

Kr 1 = <r 2 + < 2 =n 2 +1 - 6 , (6.7) 

and we can use the equation (4.20) from Fendley and Intriligator's later paper J7| to check 
that this is indeed an identity. 

VII. CONCLUSIONS 

In this paper we have studied the thermodynamic limit of a family of perturbed non- 
unitary superconformal field theories with N = l and central charge c = — 3(4n + 3)/(2n + 2). 
We identified the two-body S-matrices of those theories with the Boltzmann weights of an 
eight-vertex model and we have seen that supersymmetry implies the free fermion condition 
and the value H = — 1 for the magnetic field of the associated XY model. This allowed us 
to diagonalize the transfer matrix by Felderhof 's method and to develop its TBA. We found 
the algebraic equations that the TBA leads to in the UV limit. Solving those equations led 
to the correct value for the central charges. This analysis has thus confirmed the validity of 
the S'-matrices ( p.3|) ( p75|) and, in particular, ruled out possible "CDD" correction factors 
for those S'-matrices. 

Another interesting result was the verification that the N = 1 TBA systems obtained 
here can be obtained by folding in half a series of N = 2 supersymmetric TBA systems. This 
idea was first proposed by Melzer in |J. 
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There are several possibilities for further work. We have seen that Melzer's proposal 
works and that the folding idea can indeed be extended to TBA systems for non-diagonal 
S'-matrices. It should be possible to understand this folding directly at the level ^-matrices, 
possibly as some generalization of the non-critical orbifold construction of P. Fendley and 
P. Ginsparg [Q. This will be subject of a future publication [|TT|]. 



Recently a lot of (deserved) interest has been given to 1 + 1 field theories with a boundary, 
i.e., defined on a half-line. It is not very difficult to find a boundary S-matrix or reflection 
matrix R(9) for the theories that we discussed in this paper. In particular, we can easily 
find the boundary analogue of the Bose- Fermi .S-matrix S BF (9), see ( |2.5|) . In order to do 
that we have to solve the equation 

Q B (9)R BF (9) = R BF (9)Q B (-9) , (7.1) 

where Q B {9) is the supersymmetry charge on the half-line. In terms of the bulk supersym- 
metric charges Q(9) and Q{9) g| we have Q B {9) = Q(9) + Q(9). This gives 

Q B (0) = e*(T x + e-$(T v . (7.2) 



Solving the condition (|7.1j) we obtain 



R BF {9) = Z{9) 



I cosh(f + %\ ) » 



(7.3) 



y cosh(f - %\ ) j 

where Z(9) is a normalization factor. We have verified that this boundary matrix R BF (9), to- 
gether with the bulk S'-matrix S BF (9), satisfies the boundary Yang-Baxter equation. Clearly, 
we can extend the entire analysis performed in || and in this paper to the boundary case 
and explore the physical meaning of such results. This will be the subject of future work. 
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